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The representation of quark distribution functions in terms of nonlocal operators is combined with 
a simple diquark spectator model. This allows us to estimate these functions for the nucleon ensuring 
^Q ' correct crossing and support properties. 

<D " 

(^ I I. INTRODUCTION 



Quark distribution functions appear in the field-theoretical description of hard scattering processes as the parts 
that connect the quark and gluon lines to hadrons in the initial state. These parts are defined through (connected) 
qq ', matrix elements of nonlocal operators built from quark and gluon fields. The simplest, but most important ones, are 
£SJ ■ the quark-quark correlation functions p|-0J . For each quark flavor one can write down 

<n: 4 

^ $ ij (p,P,S) = J-^±e- i ^ (P,S\^(OM0)\P,S), (1) 

The hadron states are characterized by the momentum and spin vectors, \P, S). The quark momentum is denoted p. 
O ! A summation over colors is understood in <&. 

In a particular hard process only certain Dirac projections of the correlation functions appear. For inclusive lepton- 
hadron scattering, where the hard scale Q is set by the spacelike momentum transfer — q 2 = Q 2 , the correlation 
function $ appears in the leading order result in an expansion in 1/Q and the structure functions can be expressed 
r* ' in quark distribution functions. 

The hard momentum scale q and the hadron momentum P define the lightlike directions n±. In inclusive lepton- 
S^ i hadron scattering they are related to the hadron momentum P and the hard momentum q as 

•^ P = ^n„ + Q n+ , (2) 

q = T2 n -~T2 n+ - (3) 

The quark momentum p and the spin vector S are also expanded in the lightlike vectors, 

xQ x b {p 2 +p 2 T ) 

P= ^V=2 n++ X QV2 n ~ +PT > (4) 

S = ;= Tl + -=r- Tl- + 5t- (5) 

MxbV2 QV2 

Thus x represents the fraction of the momentum in the + direction carried by the quark inside the hadron. The 
lightlike components a = a ■ n^ are defined with the help of the vectors n± . 

It should be noted that there are other contributions in the leading cross section resulting from soft parts with A + 
gluon legs. These can be absorbed into the correlation function providing the link operator that renders the definition 
in Eq. |l| color gauge invariant. Choosing the A + = gauge in the study of $ the link operator reduces to unity. 
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Although hard cross sections can be expressed in terms of distribution and fragmentation functions, these objects 
cannot be calculated from QCD because they involve the hadronic bound states. We follow here a different route. We 
investigate the structure of lightcone correlation functions (including the effects of transverse separation of the quark 
fields), we illustrate the consequences of various constraints and we estimate the distribution functions in a simple 
model. Particularly suitable is a model in which the spectrum of intermediate states that can be inserted in Eq. |l| is 
replaced by one state, which if the hadron is a baryon will simply be referred to as a diquark. 

II. QUARK CORRELATION FUNCTIONS 

The most general expression for $ is j|,[| : 

$(p, P,S) = MA X + A 2 f + A 3 i> + MA 6 ^ 75 + ^p ■ S P 75 



+^p-S^ l5 +iA 9 a" u ^S ti P u 



M 
+i A 1Q 0^ l5 S^ + i^ P -S <j» v 75 p M P v , (6) 

where the amplitudes Ai depend on a = 2p ■ P and t = p 2 . Hermiticity requires all amplitudes Ai(<r, r) to be real. 
Time reversal invariance has also been used. 

As the variation along n_ up to order M 2 /Q 2 is irrelevant in a hard process, one encounters the quantities 



$[ r l(x,p T )= l -J dp-^r^Y) 



(7) 

p+=xP+ 



The projections of $ on different Dirac structures define distribution functions. They are related to integrals over 
linear combinations of the amplitudes. The projections which contribute at leading order are 

^ + ^(x,p T ) = \g lL (x,p 2 T ) + E2j^g 1T (x,p 2 T ), 

<+ ~" 1 ' ^ -" - ' ,2\ , Pt f^ h x r i\ , Pt'St ,!_ 



^ Tl&pr) = S l T h 1T (x,p 2 T ) + jL \\ht L {x,p 2 T ) + ^^ht T {x,p 2 T ) 

After integration over the quark transverse momentum, these leading twist functions have well known probabilistic 
interpretations. The function fi(x) gives the probability of finding a quark with light-cone momentum fraction x in 
the + direction (and any transverse momentum). g\(x) is a chirality distribution: in a hadron that is in a positive 
helicity eigenstate (A = 1), it measures the probability of finding a right-handed quark with light-cone momentum 
fraction x minus the the probability of finding a left-handed quark with the same light-cone momentum fraction (and 
any transverse momentum). h\(x) is a transverse spin distribution: in a transversely polarized hadron, it measures 
the probability of finding quarks with light-cone momentum fraction x polarized along the direction of the polarization 
of the hadron minus the probability of finding quarks with the same light-cone momentum fraction polarized along 
the direction opposite to the polarization of the hadron. 

The subleading projections can also be expressed in distributions. The integrated distributions are e(x), grix) and 
hz,(x). Important relations are 

9t(x) = gi(x) + — g$(x), (8) 

h L (x) = h 1 {x)-£hH 1) (x)„ (9) 

where 

d 2 PT ygjj f{ X ,p%). (10) 

The functions 32 = <?t — ffi and /12 = 2(/il — hi) thus satisfy the sum rules 



1 

dx g 2 (x) = -<7$(0), (11) 

o 
l 

dx h 2 {x) = 2/i^ 1) (0), (12) 



If the functions <^ T and /i^ vanish at the origin, we rediscover the Burkhardt-Cottingham sum rule H and the 
Burkardt sum rule 0. 

III. THE DIQUARK SPECTATOR MODEL 

The basic idea of the spectator model is to treat the intermediate states that can be inserted in $ in Eq. [j] as 
a state with a definite mass Mr, called the spectator. In other words, we make a specific ansatz for the spectral 
decomposition of the correlation function. The quantum numbers of the intermediate state are those determined by 
the action of the quark field on the state \P, S 1 ), hence the name diquark spectator. The inclusion of antiquark and 
gluon distributions requires a more complex spectral decomposition of intermediate states. We restrict ourselves to 
the simplest case. $ is then given by 

*3 = 7A3 {P,S$j(0)\X W ) 8(Pr) 8 [iP-Pf Mr] (X^\M0)\PS) (13) 

where Pr = P — p and X^ represents the spectator, whose possible spin states are indicated with A. 

The spin of the diquark state can be (scalar diquark s) or 1 (axial vector diquark a). The matrix element appearing 
in the RHS of (Oh is given by 

(X S MQ)\P,S) = f^ 1 -) Tfc Ut(P,S), (14) 

\p-mj lk 

(x^\Mo)\PS) = e;^(j^-y r$ut(p,s), (15) 

for a scalar and a vector diquark, respectively. The matrix elements consist of a nucleon-quark-diquark vertex T yet 
to be specified, the Dirac spinor for the nucleon Ui(P,S), a quark propagator for the untruncated quark line (m is 
the mass of the quark) and a polarization vector e^ in the case of an axial vector diquark. The next step is to fix 
the Dirac structure of the vertex T. We assume that 



? s = 1 9s(p 2 ), (16) 

9a(p 2 ) ( M , P" 



T° = ^ 75 (7" + ^7), (17) 



where gnip 2 ) (with R being s or a) are form factors that take into account the composite structure of the nucleon 
and the diquark. In the choice of vertices, the factors and projection operators are chosen to assure that in the target 
restframe the vector diquark states are pure spatial and the nucleon spinors have only upper components, in which 
case e.g. the axial vector diquark vertex reduces to X N cr ■ e\ q - With our choices we obtain 



I 2 ' 
* K = ^rrf — ,Y "\,2" BJ ii> + m)(r + M) (1 + a Rl5 $) tf + m), (18) 



\ g R (p 2 )\ 2 S(t-o- + M 2 - M 2 r , 

2(2^) 3 (p2 _ m 2)2 



where a s = 1 and a a = —1/3. In obtainig this result we used as the polarization sum for the axial vector diquark 

Yl\ e ^ € » = — 9v + P^Pu/M 2 , which is consistent with the choice that the axial vector diquark spin states are 
purely spatial in the nucleon restframe. We will use the same form factors for scalar and axial vector diquark being 
of the form 

2 
T — Tft 

A is another parameter of the model and N is a normalization constant. This form factor has the advantage of killing 
the pole of the quark propagator M . 



IV. RESULTS AND DISCUSSION 



Using Eq. Hq we can compute the amplitudes Ai shown in Eq. pi Taking out some common factors by defining 

N 2 



At 



2(2vr) 



5 (t - a + M 2 - M R ) 

|t-A 2 | 2q 



A> 



(20) 



we obtain A4 — A$ = A12 = and 



A, = -A 6 /a R = g ((M + to) 2 - Af 2 ) + (r - to 2 ) (l + £) , 

A 2 = -Ag/a R = - (r - to 2 ) , 

is = -Aw/ait = (Af + to) 2 - Af| + (r - to 2 ), 

A 7 = 2 ai? mM , 

I 8 = -In/aR = 2a il M 2 . 



Introducing the function 



A|(») = A 2 (l -x) + xM R - x(l - x)M 2 



one gets, e.g., for fi(x,p 2 -,) the result 



fi(x,Pr) 



N 2 (l-x) 2a - 1 (xM + mY+p 2 T 



16tt 3 



(p 2 T + X 2 R f a 



(21) 



(22) 



Although there is a certain freedom in the choice of the parameters, one immediately sees that the occurence of 
singularities in the integration region will cause problems that can be avoided if we require A^(a;) to be positive, what 
imposes the condition M R > M — A. Provided this condition is fulfilled one obtains the integrated distributions: 



JV 2 (1 - a:) 2 "- 1 
h[X) 32tt 2 (a - l)(2a - 1) 



2(a - I) (xM + m) 2 + \ 2 R {x) 



2a-l 



\2a-l 



N 2 a R (l-x) 2 
9l[X) 327r 2 (a-l)(2a-l) 



(AIW) 
2(a-l)(xM + TO) 2 -A| e (a;) 



/ii(x) = 



_ Ar 2 ai? , (1 - a) 



2q-1 



16^ 2 (2a-l) (A 2 fl (x)) 



(Mi + to) 2 

2q-1 



Explicit expressions for the twist three functions can also be written 
An example of a p|,/2M 2 -weighted distribution is 



,(!) 



9it ( x ) 



N 2 a R (1 - x) 



2a- 1 



(xM + m)M 



32^ 2 (a-l)(2a-l) {Xr (x)) 



2a-2 



(23) 



Up to now, we have not specified flavor in the distributions. For the nucleon we distinguish two types of distributions, 
/-f and fi , etc. This corresponds to the spectator having spin or spin 1, respectively. In that case spin diquarks 
are in a flavor singlet state (S) and spin 1 diquarks are in a flavor triplet state (T) in order to combine to a symmetric 
spin-flavour wave function. Since the coupling of the spin has already been included in the vertices, we need the flavor 
coupling 



^ = 71 ' u s °^ + 7l' u T °^ ~ ~/3 ' d Tl ^' 



(24) 



to find that the flavor distributions are /" = | f* 



\ /f) ft — fi an( i similarly for g\, hi, e, gr and h^. The 
proportionality of the numbers is obtained from Eq. |24], while the overall factor is chosen to reproduce the sum rules 
for the number of up and down quarks if /* and /" are normalized to unity upon integration over p T and x. This 
will fix the normalization N in the form factor. The factors a s — 1 and a a — — 1/3 in the distribution functions will 
produce different u and d weighting for unpolarized and polarized distributions. Further differences between u and d 



distributions can also be induced by different choices of Mr, A or a. We take for the nucleon a = 2 to reproduce the 
right large x behavior of /", i.e. (1 — x) 3 , as predicted by the Drell-Yan-West relation and reasonably well confirmed 
by data. We refrain from tuning the large x behaviour of ff to match the (1 — x) 4 form indicated by data. Since 
fi is only affected by vector diquarks, this could be easily obtained by choosing a different form factor for the latter. 
We feel that this kind of fine-tuning would take things too far with the simple model we use. Similarly, we will 
only consider one value of A. We will, however, allow for different masses for scalar and vector diquark spectators. 
The color magnetic hyperfine interaction, responsible for the nucleon-delta mass difference, will also produce a mass 
difference between singlet and triplet diquark states. Neglecting dynamical effects, group-theoretical factors lead to a 
difference M a - M s = 200 MeV. 

Another important constraint comes from the axial charge of the nucleon 



g A = dx [g?(x) - g((x)] = dx 
Jo Jo 



2 9iW + 2 Sito 



(25) 



The sensitivity to the parameters Mr and A is best illustrated by considering some characteristic values. We take 
a quark mass of 0.36 GeV, two different values for Mr (0.6 and 0.8 GeV) and three values for A (0.4, 0.5 and 0.6 
GeV.) In Table B the values of some moments are given. 



TABLE I. The second moment of /i, (x) = J dxxfi(x) and the first moments g\ 
given for two diquark masses and for three values of A. 



J dxgi(x) and hi = J dxhi(x) 



A (GeV) 


(X)* 


Mr = 0.6 GeV 
9? 


hf 


{x) R 


Mr = 0.8 GeV 
3? 


h? 


0.4 
0.5 
0.6 


0.367 
0.373 
0.383 


0.922 a R 
0.792 a R 
0.664 a R 


0.965 a R 
0.896 a R 
0.832 a R 


0.232 
0.255 
0.278 


0.643 a R 
0.524 aR 
0.412 a R 


0.821 a R 
0.760 a R 
0.707 a R 



Fig. [I] shows the twist two distributions for different values of the mass of the spectator and of the parameter A. 
We can see that an increase of Mr induces a shift on the peak of the valence distribution f\(x) towards lower values 
of x and a decrease in its second moment. This is due to the fact that a spectator with higher mass corresponds to 
more massive intermediate states which can contain sea quarks. In these circumstances the momentum carried by 
valence quarks decreases. In the non relativistic limit the equality f\{x) = gi(x) = h\{x) holds, a result that could 
be anticipated. 

We now use the axial charge of the nucleon to find the most suitable values for A and the diquark masses. The 
values A = 0.5 GeV, M s = 0.6 GeV and M a = 0.8 GeV give gA = 1.25, close to the experimental value. Fig. shows 
the distribution f\{x) multiplied by x for these values of the parameters. We find a satisfactory qualitative agreement 
with the valence distributions of Gliick, Reya and Vogt [[10| . 

Finally, Fig. shows the distributions g%(x) and g^ix), for which we find a small violation of the Burkhardt- 
Cottingham sum rule, in agreement with Eq. [Ll[ 
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FIG. 1. Twist two distributions for the nucleon. The plots on the top represent fi(x), the ones on the middle show gi(x)/an 
and at the bottom we have h\(x)/an. The plots on the left correspond to A = 0.4 GeV and the ones on the right to A = 0.6 
GeV. The full line corresponds to M R = 0.6 GeV and the dashed line to Mr = 0.8 GeV. 
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FIG. 2. Twist two distributions for the nucleon. The plot on the left shows xf^(x) (full line) and xff{x) (dashed line) for 
Ms = 0.6 GeV, M a = 0.8 GeV and A = 0.5 GeV. The plot on the right shows the "low scale" valence distributions of GRV. 




FIG. 3. Distributions gl (x) (solid line) and g$(x) (dashed line) for M s = 0.6 GeV, M a = 0.8 GeV and A = 0.5 GeV 



